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Abstract 

Virial relations for the Dirac equation in a central field and their applications to 
calculations of H-like atoms are considered. It is demonstrated that using these re- 
lations allows one to evaluate various average values for a hydrogenlike atom. The 
corresponding relations for non-diagonal matrix elements provide an effective method 
for analytical evaluations of infinite sums that occur in calculations based on using 
the reduced Coulomb-Green function. In particular, this method can be used for cal- 
culations of higher-order corrections to the hyperfine splitting and to the g factor in 
hydrogenlike atoms. 

1 Introduction 

In non-relativistic quantum mechanics, the virial theorem for a particle moving in a central 
field V{r) is given by the well known equation: 

where (T) denotes the average value of the kinetic energy in a stationary state. This theorem 
can easily be derived from the equation [1] 

|(A)4([i/,A]).0, (2) 

if we take A — (r-p). Equation (2), which is generally called as the hypervirial theorem, can 
also be employed to derive virial relations for diagonal matrix elements of other operators [2] . 
An extention of these relations to the case of non-diagonal matrix elements was considered 
in [3, 4]. 

For the Dirac equation in a central field, the virial theorem was first derived by Fock [5]. 
If we denote the bound-state energy by E, it gives 

E = (mc^P) + (r^> + (V) , (3) 

where (3 is the Dirac matrix. For the Coulomb field, one easily finds 

E = {mc^f3) . (4) 
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Some additional virial relations for the Dirac equation were obtained by a number of authors 
(see, e.g., [2, 6] and references therein). Virial relations, which yield recurrence formulas for 
various average values, were obtained in [2, 7, 8, 9]. In the case of the Coulomb field, these 
relations can be employed to derive explicit formulas for the average values (r"), (r"/?), and 
(ir*(a-n)/3), where a is a vector incorporating the Dirac matrices, n = r/r, and s is integer. 
The corresponding recurrence relations for non-diagonal matrix elements were derived in [9] . 
In the case of the Coulomb field, it was found that these relations can be employed to derive 
explicit formulas for the first-order corrections to the Dirac wave function due to interaction 
with perturbativc potentials of the form ^ r*, r*/3, r''{a ■ n), and ir^{a ■ n)/3. Later on 
(see references below), this method was used for calculations of various corrections to the 
energy levels, to the hyperfine structure splitting, and to the bound-electron g factor. In 
constract to direct analytical and numerical calculations, the virial relation method allows 
one to derive formulas for various physical quantities by simple algebraic transformations. 

In the present paper, following mainly to [9] , we derive the virial relations for the Dirac 
equation and examine their applications to calculations of H-like atoms. Relativistic units 
{h = c = 1) are used in the paper. 



2 Derivation of the virial relations for the Dirac equa- 
tion 

For the case of a central field V(r), the Dirac equation has the form 

{-ia -V + f3m + y(r))V'(r) = E^P{r) . (5) 
The wave function is conveniently represented by 

^{r)=( , (6) 

V «/('')"-Km(n) / ' ' 

where k = (— l)-'^'^^^'^(j + 1/2) is the quantum number determined by the angular momen- 
tum and the parity of the state. Substituting this expression into (5), we obtain the radial 
Dirac equations 

^ + -G-{E + m-V)F = 0, (7) 
dr r 

^-lF+{E~m~V)G = 0, (8) 
dr r 

where G{r) = rg{r) and F{r) = rf{r). Introducing the operator [10] 

H^ = -iay^ + a^- + (J^m + V , (9) 

dr r 

where Ux, (Jy, and are the Pauli matrices, and denoting 

_ / G(r) 



<t>ir) = ( ) , (10) 
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we obtain 

H,<t> = E4>. (11) 

The operator is self-adjoint in the space of two-component functions satisiying the bound- 
ary conditions 

(/.(O) = 0M=O. (12) 
The scalar product in this space is defined by 

/■oo 

{a\b)= / dr{GaGb + FaFb). (13) 
Jo 

Let us denote the eigenvalues and the eigenvectors of the operator by En^, and (f)nKi 
respectively, where n is the principal quantum number. Taking into account the self- 
adjointness of iJ^, we can write down the following equations 

{n'K'\{H,,Q-QH,)\nK) = {E^, - E^,){n' n'\Q\nn) , (14) 
{n'K'\{H^,Q + QH^)\nn) = {En> ^> + En^){n' n'\Q\nK) . (15) 

Substituting Q = r*, icr^r* into equation (14) and Q = Ozt'^ ■, Oxt" into equation (15), and 
using the commutation properties of the Pauli matrices, we obtain [9] 

{En'K,' - EnK.){n'K'\r''\nK) = -s{n' K'\iayr''~'^\nK) 

+{k' - K){n'K'\axr^-'^ |n/t) , (16) 
{En'K' - Eni^){n' K'\iayr''\nK) = s{n' K'\r''~'^\nK) 

-(«' + K){n'K'\azr''-'^\nK) + 2m(n' n'la^r^lnK) , (17) 
{En'K' + EnK){n'K'\(7zr''\nK) = s{n' K'laxr^^^lnK) 

— {k' — k) {n' k' \iayr^^^ Inn) + 2m{n' K'\r'^\nK) 

+2{n'K'\azVr''\nK} , (18) 
{En'K' + EnK){n'K'\axr^\nK) = -s{n'K'\a;^r^~'^\nK) 

+{k' + K){n'K'\r'-'^\nK) + 2{n'K'\axVr'\nK) . (19) 

In the next sections, we apply these equations for calculations of the average values of various 
physical quantities as well as for calculations of various higher-order corrections. 

3 Application of the virial relations for evaluation of the 
average values 

Let consider equations (16)-(19) for the Coulomb field {V{r) = —aZ/r) and for diagonal 
matrix elements {n'n' = nn). Denoting 

/•oo 

= / drr'{Gl,+Fl,), (20) 
Jo 
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f-OO 

= / dTr^{Gl^-Fl,), (21) 

POO 

= 2 dr r'Gn.F^. , (22) 
Jo 

we obtain [2, 7, 9] 

2mA'-2En^B' = 2aZB'-^ ^ sC'^ (23) 
2mC' = -sA'-^ +2kB''-\ (24) 
2£;„«;C^ = 2kA''-i - sB^-i - 2aZC"'-i . (25) 

From these equations, one easily finds 

-[{s + 1)S„« + 2Km]A'' + [{s + l)m + 2k£;„«]jB^ + 2aZmC^ = (26) 

Using equations (23)-(26) for s = 0, 1, we obtain 

m ' aZ m? ' to^ 

In addition, according to the HeUmann-Feynman theorem, we have 

— = (riK — — InK) = {nK\ — riK) . (28) 

OK OK r 

It yields 

1 _ dEn^ _ {aZfKm 



where N = + rirY + {ot^Y^ 7 = -v/k^ — (a^)^, and rir = n — |«;|. The derivative with 
respect to k in equation (29) must be taken at a fixed Ur- Using the formulas for C", 
B~^, and given above and reccurence equations (23)-(26), we can calculate the integrals 
A^, B*, and for any integer s. Explicit formulas for these calculations were derived in 
[7]. The formulas expressing the integrals 5**+^, and C*"*"^ in terms of the integrals 

A^ B^ and C are 

= {2aZEn^{s + l)A' + 2aZm{s + 2)5" 

-(s + l)[sm + 2(m + KEn^)]C'}{2{s + 2){m^ - EIJ}~^ , (30) 
B'+^ = {2aZm{s + 1)A' + 2aZEn^{s + 2)B' 

-{s + l)[2Km +{s + 2)Er,^]C'}{2{s + 2){m'' - J}-i , (31) 

C"+i = —[2kB'-Is + 1)A']. (32) 
2m 

The reversed formulas are 

_ 4aZ(5 + 2)(toB«+i - En^A'+^) {s + l)m + 2kE^^ 



(s + l)[(s + l)2-472] aZm[(s + 1)2 - 472] 
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X {[(s + 2)m + 2KEn4B'+^ - [{s + 2)En^ + 2Km]A'+^} , (33) 
4aZ{mA'+'^ - E„^B'+^) (s + + 2/tm 



(s + 1)2 - 472 aZm[{s + 1)2 - 472] 

x{[(s + 2)En^ + 2Km]A'+^ - [{s + 2)m + 2KEr,^]B'+^} , (34) 

C = ^^^[{s + l)En^ + 2Km]A' -^^^[{s + l)m + 2KEn^]B' . (35) 
Employing these formulas, one easily finds 

^ - 2m2 ' ^^^> 

^ = :^('^ +"-^)' (37) 

rr)2 - R2 

= ^ (38) 

2 _ 2(aZ)2^[2«(7 + n,)-iV]m2 

7V4(472- 1)7 ' ^"^^^ 

0-2 2(aZ)2[27^Ar-/,(7 + n.)]m2 

^ = A^^Th ' 

2 _ 2{aZf[2^{j + nr)-N]m^ 
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2{aZfm^ 



X [7V2(i + 27^) - 3KAf(7 + n^) + 3(aZ)2(Ar2 - ^2)] . (42) 



It should be noted here that the integral A~^ exists only for \k\ > 2. The integral 
occurs in calculations of the bound-electron g factor. The integrals C^^ and A~^ occur in 
calculations of the magnetic dipole and electric quadrupole hyperfine splitting, respectively, 
(see [11] for details). Formulas (27)-(29) and (36)-(41) were also employed in calculations of 
the recoil corrections to the atomic energy levels [12, 13]. 



4 Application of the virial relations for calculations of 
higher-order corrections 

In calculations of higher-order corrections to various physical quantities one needs to evaluate 
the sums 

, , ''^'^'sr~^"''^\n'K,')(n'K'\R^\nK) ,,„, 

\t,S,K',nK)= Yl F _F ' (43) 

n' 

where i?f = r*, i?2 = f^zT^ , -R3 = (^xi"^ 1 and R% = iayV^ . For instance, to derive the first-order 
correction to the hydrogenic wave function due to the magnetic dipole hyperfine interaction. 
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we need to evaluate the expression (43) for R^^. Let us consider how virial relations (16)- 
(19) can be employed for calculations of these sums in the case of a hydrogenlike atom 
{V{r) = —aZ/r). For this case, equations (16)-(19) can be rewritten in the following form 

(£„„-i;„,«,)^^'K',n« = -'2rnC^'.',n.-sAl7^',n. + i^' + '^)K~.\n.^ (45) 

+(«' - n)D:^7Xn. + 2«^S:r;, _ , (46) 
+2«^t^:r;,„,, (47) 

where 

/•OO 

), (48) 

^0 

K'.',n. = / drr^{Gn'.'Gn.-Fn'.'Fn.), (49) 
Jo 

/■oo 

C^'.',n. = / drr'iGn'.'Fn^ + Fn'.'Gnn), (50) 

Jo 

/>oo 

K'.',n. = / drr%Gn'.'Fn,-Fr,,,,Gn.). (51) 
Jo 

Prom equations (44)- (47), we obtain 

{EriK - En'K'){EniiDn'K',nK + '^^n' k' ,711^) 

= [sEn^ - m{K' + k)]<7,\„, + [sm + E^^in' + n)]B'J^,^n. 
+2aZmC:j-Xn,. (52) 

For riK = u'k', this equation turns into equation (26). 

Let us consider first the Ccise K — /^ • Taking into account that ^°/^ „^ = Sn'n, we obtain 

En'K,nK ~ ~iEnK ~ En'K){EnKD„if^ nK ''^^n'K,nK 

+«^^°'«,n. - «S°'«,nJ + ^^nW • (53) 

Multiplying this equation with |n'«;) and summing over n', we derive 

|2,0,«;,n«;) = —{I — \nK){nK\) 
m 

x{Enf^i(jyr + mcTxr + aZiuy — Kaz)\nK) , (54) 

where I is the identity operator. From equations (52) and (54), we can derive the sum 
|3, 0, K, UK). Then, using equations (45)-(47) and (52), we can calculate all the sums \i,s, K,nK) 
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for i = 1,2,3 and s = 0, 1, 2, ... . In particular, for the sum |3, 1, k, uk) that occurs in calcu- 
lations of various corrections to the bound-electron g factor, we find 



|3, 1, K, UK) 



K 

m 

X 



2(1- \nK){nK\) 

(TTl \ 



\nK) . 



(55) 



The sums |4, s, k, n«;) for s 7^ —1 can be calculated by employing equation (44). The sum 
|4, —1, K,nK) is easily derived from the relation 



{Enii - En'K,){n'K\lnr\nK) , 



which is obtained by differentiation of equation (44) with respect to s. 
Let us consider now the case k' ^ k. From equation (44) we find 



|3, —1, k', nn) = 



K — K 



j\nK). 



(56) 



(57) 



Then, using equations (44)-(47) and (52), we can calculate all the sums \i, s, k' ,nK) for 
i = 1,2,3,4 and s = —2,-3,-4,... (if, of course, the corresponding sum exists). For the 
sum |3, —2, k', tik) that occurs in calculations of various corrections to the hyperfine splitting, 
we find 



|3, —2, k', tik) 



x[[l-(- + «'f](|^ + (-'--)--^ 



K)m — En 



+ - 



K — K' 



For K = ±k', the corresponding sums can be calculated by taking the limit k' 
taking the limit k' — > «; in equation (44), we obtain 

3, —1, K, tik} = — tik) , 

OK 



(58) 

±K. So, 

(59) 



where, as in (29), the derivative with respect to k must be taken at a fixed n^. Then, 
the other sums with k' = k and s = -1,-2,-3,... can be calculated by using equations 
(44)-(47) and (52). In particular, we obtain 



\3,—2,K,nK) = 



^{aZf + (1 - 4k2) 

2{aZ)^Km 



2aZaxT ^ + ^.aZKiGyr ^ 



(1 — 4^2) 



\nK) 



K 

2aZ{2E„ 



'V 



A{aZy + (1 - 4k2) Ok 



tik) . 



(60) 
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The case k! = —n can be considered in the same way. 

Concluding this section, we give the expUcit formulas for -^\nK) for the Is and 2s states 
[14]. For the Is state: 



^1 \ 



GUr) 
Fisir) 



where 



Gu{r) = 



Fu{r) 



1 t 1 

k 



VT+7 



exp(-i/2)F( 
exp ' 



V;(27+l) 



+ (7 + 1) 



(-i/2)t^( 



^(27+1) 



+ (7+1) 



1 t 1 , \ 
+ 2:^-2-7'^')' 

2aZmr , (2aZ)im^ 
t = — — — , fc — 



N ' 2v/2r(27+l) ' 

r(a;) is the gamma-function, and il>{x) = ^ lnr(a:). For the 2s state: 



^1 \ 



G2s{r) 



where 



G2s{r) = 



\2(N - 



2' iV2 - 2 l2(iV- 1) 
'2N'^ -4N^ + 5N'^ -3N + 2 2^^(27+1) 



+ 



+ 



V 2N{N - 1)2 

N^~2N^ + N~2 

+ 27VV'(27 + 1) 



N 



N-1 



2{N -1) 
-tlnt-2Nlnt 



F2s{r) = 



}■ 

l2fiV- 



-^'^"P(-^)*'^-fl2(iV-l) 

/2iV4-2iV3 + iV2 + 3iV-2 ^ 2^^(27 + 1) 
~l 2iV(iV- 1)2 

+ 57V2 - 87V - 4 



+ 



TV 



+- 



2N{N-1) 



+ 2(7V + 2)V(27 + 1) 



+-^^— Yilnt-2(7V + 2)lni} , 



(61) 



(62) 



(63) 



(64) 



(65) 



(66) 



8 



2aZmr 
N 



k' 



r(27 + l) ^8N{N + l)\ N 



^2aZm 



5 Calculations of the bound-electron g factor and the 
hyperfine splitting in H-like atoms 

For the last few years a significant progress was achieved in calculations of the bound-electron 
g factor and the hyperfine splitting in H-like atoms. Formulas (55) and (60) were extensively 
employed in these calculations. 

In [15], a complete a.Z-dependcncc formula for the recoil correction of order m/M to the 
g factor of an H-like atom was derived. According to this formula, which was confirmed by 
an independent derivation in [16], the recoil correction is given by (e < 0) 



A5 



1 



duj 



_d_ 



{a\\p'^-D\uj) + eA'',,] 



xG{uj + E^)\p^-D\uj) + eA'l,]\a) 



(67) 



H=0 



Here jjbQ is the Bohr magneton, rua is the angular momentum projection of the state a, 
p'^ = — iV*^ is the momentum operator. Ad = [Ti x r]/2 is the vector potential of the 
homogeneous magnetic field Ti, directed along the z axis, D^{uj) = —4:'KaZa^D''''{u), 



47r L r 



'i\uj\r) 



2 



(68) 



is the transverse part of the photon propagator in the Coulomb gauge. The tilde sign 

indicates that the related quantity (the wave function, the energy, and the Coulomb-Green 
function G{io) = Ylin ~ -^n(l ~ *0)]~^) must be calculated at the presence of the 

homogeneous magnetic field H directed along the z axis. In equation (67) and below, the 
summation over the repeated indices (fc = 1, 2, 3), which enumerate components of the three- 
dimensional vectors, is implicit. For the practical calculations, this expression is conveniently 
represented by the sum of the lower-order term and the higher-order term, H^g = AgiL + Ag'H, 
where 



A<7L 



1 1 

HOTUa 2M 

1 m 



A(a|{p2 _ ^[(a . p) + (a . „)(„ . p)]}la) 



J H=0 



(69) 



cG(c. + ^„)(D'^(a;) + ^^)la) 



(70) 



H=0 
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where V{r) = —aZ/r is the Coulomb potential induced by the nucleus and n = r/r. 
The lower-order term can be calculated analytically by employing formula (55) and the 

formulas for the average values presented above. Let us consider this calculation in details 
[15]. According to equation (69), the lower-order term is the sum of two contributions, 



AgL = ^gj^^ + ^q'l'- The first contribution is 



,(2) 



Ag. 



1 



1 



{6a\ 



— {a • p + (a • n)(n • p) la) 
r J 



(71) 



where Sa is the first-order correction to the electron wave function due to interaction with the 
homogeneous magnetic field. Taking into account that p^ = (a-p)^ and (a-p) = H—fim—V, 
one easily obtains 

(5a|p2|a) = {5a\{E^+(im-V){E^-(im-V)\a) 

+i{5a\{cx-VV)\a) . (72) 

The second term in equation (71) can be transformed as (see, e.g., [13]) 



cxZ 

-{5a\ [a. • p -I- (a • n)(n • p)]|a) 

r 



-{5a 



aZ 



2Ea - 2(3m - 2V 



+- (a-n)(/3«; + l) \a) . 
r J 



The wave function correction \5a) is defined by 



\6a) 



E 



Ea — En 



(73) 



(74) 



Since the operator sandwiched between \a) and \6a) in equation for Ag^"^ conserves the 
angular quantum numbers, we need only that component of \Sa) which has the same angular 
quantum numbers as the unperturbed state \a). Using formula (55), we easily find 



\6a)^ 



ir(r)l^_KT»a(n) 



where 



X{r) 

Y{r) 
bo 



- m 

- M[ 



2mK — m + 2K£Ja ol'Z' 
T— 2 r^ 

2mK + m — 2At£Jo aZ i 
— r 



2to2 



5(0 + 



K- 2«;2 

K + 2k^ 
2m? 



<7(r)}, 
f{r)}, 



K 



2 i(i + l)' 



(75) 

(76) 

(77) 
(78) 



g{r) and f{r) are the radial parts of the unperturbed wave function defined above. Inte- 
grating over the angular variables in equations (72) and (73), we find 

^^L^ = j^^^fdrr'{x{r)g{r)[-2Vm-V' + El-m'] 
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+Y{r)f{r)[2Vm - + - 
+ [X(r)/(r)+F(rMr)]^Ac}. 

Substituting expressions (76) and (77) into equation (79), we obtain 

A^W = _^^/aZ?^^^^C° + (aZ)2^C7-^ 
i(j + l)Ml m ^ ' m 



(79) 



3mK-2K^Ea^. 
2^^ ^ 



-1 



(80) 



Using the explicit expressions for the integrals A^, , and C* given above, we obtain 

(1) ^ m - 

M 2j{j + 1) m2 ■ 

Consider now the contribution Ag^ ' : 

=-;;^M<'^lO^-^I^^"l^)|'^>• 



(8l) 



(82) 



Integrating over the angular variables and employing the explicit expressions for the radial 
integrals derived above, we obtain 



^9^' = - 



m 
M 2j{j 



m + Ea 



+ (2j-0(2j-/ + l)^^^:^-«2 



2m 

m2 - E? 



2m 



For the sum A^il = ■^^l^^ + we find 



m 2k^E^ + nmEa — 



M 2m2j(j + l) 
To the two lowest orders in aZ, we have 



A5L 



m 



M j{j + 1) 
For the Is state, formula (84) yields 



K\ {aZf 



. m , ^.n m 
Ag^ = -iaZr- 



{aZY 



Mz[l + ^l-{aZff 



(83) 



(84) 



(85) 



(86) 
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The first term in the right-hand side of this equation reproduces the result of [17, 18]. The 
higher-order term A^h was evaluated numerically for the Is state in [19]. Formula (55) was 
also extensively used in that evaluation. 

Let us consider now the derivation of the nuclear-size correction to the g factor of a 
low-Z H-like atom [20]. To find this correction, we have to evaluate the expression 

2 '^ {a\5V\n){n\\e\a-Aa\a) 

^9=J^}^ ' ^''^ 

where SV determines the deviation of the potential from the pure Coulomb one. Integrating 
over the angular variables, we obtain 

2Km ""s—^^ {nK\6V\n' K}{n' K\rax\nK) 

'^1UTT)Z-. Er.. - E^,. ' 

where \nK) are the two-component radial wave functions defined above. Substituting ex- 
pression (55) into this equation, we get 

2/^^ r r / Tfx \ "I 

^9 = . , i {nK\5V [EnK - — ) riay + mrax + aZiay - na^ \nK) 
j{j + l)m^ LV 2k/ m j I / 

— {nn\ ^(yEjiK. — ficTy + mrcFx + aZiay — Kaz \nK) 

x{nK\6V\nK)Y (89) 

We assume that the nuclear charge distribution is described by a spherically symmetric 
density p{r) = p{r), which is normalized by the equation 



dr p(r) = 1 . (90) 
The Poisson equation gives 

A{SV){r) = ATTaZ[p{v) - (S(r)], (91) 

where A is the Laplacian. When integrated with 5V, the radial functions g{r) and /(r) can 
be approximated by the lowest order term of the expansion in powers of r. It follows that 
we have to evaluate the integral 



Using the identity 



dr r^r^''-^5V . (92) 
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and integrating by parts, we find 

1 







dr r 



27(27+1) 



-I- 



Ar^T SV= dr 



1 



27(27+1) 



dr r^ r^'^ p{r) 



27(27+1) 
aZ 



27(27+1) 



r^'' A{5V) 

Y /nuc 5 



where 



(r'^)nuc = J drr^Xr). 



(94) 



(95) 



For the correction to the g factor we obtain 



Ag = 



r(27 + 1 + n,)22T-i 



i(j + 1) 7(27 + l)r2(27 + l)nr\{N - k)7V27+2 



[nl + (TV - Kf] 1 - 2k^ - 2nr{N - «) ^ - 2k 



m 



En 



m 



xiaZf-'+^m^^r^''), 



(96) 



For ns-states, which are of particular interest, the expansion of this expression to two lowest 
orders in aZ yields 



A5 



^(aZ)^m^(r^)„uc 



i + {azy\- + 



1 12n^ - n - 9 



4n2(n + 1) 



+2V'(3) - V(2 + n) - 



(r^ ln(2aZmr/n))r 



where = ^ liir(x). For the Is state, we have 



Ag=-{aZfm'{r% 



1 + {aZf 2 - C • 



(r^ ln(2aZmr))i: 



(97) 



(98) 



where C = 0.57721566490 ... is the Euler constant. In the non-relativistic limit, we find 



for ns states and 



Ag=^{aZ)^m'{r% 



Ag=^-^^^^{aZrm^r% 



(99) 



(100) 
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for npi states. In the case of the Is state, the expression (99) coincides with the related 
formula in [21]. A similar derivation for the nuclear-size correction to the hyperfne spliting 
was performed in [22]. 

In [23, 24], formulas (55) and (60) were employed to evaluate analytically the one-loop 
vacuum-polarization (VP) corrections to the bound-electron g factor and to the hyperfine 
splitting in \ow-Z H-likc atoms. These corrections arc described by diagrams presented in 
Fig. 1 (a-c) , where the dashed line ended by a circle indicates the interaction with a magnetic 
field. According to [23], for the Is state, the VP correction to the g factor calculated in the 
Uehling approximation is given by 



16, 5tt , 



^ln(2aZ)-^)(aZ)^ 



+ 



(-^ln(../2)-^)(..)^... 



9 ^ ' ' 216 
The Uehling correction to the Is hyperfine splitting is [24] 



(101) 



8 
15 



- — ln(2aZ))(aZ)2 



/5397r IStt , , „s<t 
V 288 24 ^ ' 'jy ' ^ 



(102) 



where Ep is the Fermi energy. In [25, 26], the one-loop self-energy corrections to the hyperfine 
splitting and to the bound-electron g factor (Fig. 1 (d-f)) were calculated numerically to 
a high accuracy. Formulas (55) and (60) were also employed in these calculations. In [27], 
formula (55) was used to evaluate the nuclear-polarization effect on the bound-electron g 
factor. 



6 Other applications of the virial relations 

Applications of the virial relations are not restricted only by H-like atoms and by the pure 
Coulomb field. In [14], formulas (58) and (60) were employed to calculate the interelectronic- 
intcraction corrections to the hyperfine splitting in Li-like ions. In [29], virial relations 
(16)-(19) with V{r) ^ —aZ/r were used to evaluate the recoil correction to the Lamb shift 
for an extendc;d nuc;lcus. As an example, let us demonstrate how the virial relations with 
V{r) ^ —aZ/r can be employed to evaluate the nuclear-size correction of the lowest order 
in mR {R is the nuclear charge radius) to the integral 

/•oo 

= 2 / dr r-'^GF (103) 
Jo 

for an ns state, which we denote by |ns). This integral occurs in calculations of the hyperfine 
splitting. We consider virial relations (16)-(19) for the case n = n', k = k' and for the 
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potential 



V{r) = -— + 5V{r), (104) 
r 

where 5V determines a deviation of the potential from the pure Coulomb one due to the 
finite-nuclear-size effect. Using notations (20), (21), and (22), where G and F are calculated 
for potential (104), from equations (16)-(19) we derive 

^-2 _ = "regular terms" , (105) 

C^^ — 2{ns\(j ;,r~^V\ns) = "regular terms" , (106) 
2^-2 _ _ 2{ns\a^r-^V\ns) = "regular terms" , (107) 

where by "regular terms" we denote terms which have the nuclear-size corrections of order 
(mi?)2 and higher (it means that their integrands have more regular behaviour at r ^ 
than the integrand of has). Prom these equations, we easily obtain 

= Q A( ^.^ {ns\{2,a,r-^ -2aZa,r-^)5V\ns) 
6 — 'l[aA y 

+" regular terms" . (108) 
To the lowest order in mi?, it follows 

5C-^ = 3 _ ^^^^^2 ("^1(3^^^"' - '2aZa^r-^)5V\ns) . (109) 
An analytical evaluation of this expression to the two lowest orders in aZ yields [28] 

4 



= -^{aZ)^m^r}r,^A 1 + (aZ)' 2^{S) - V'(n + 1) 
n'^ I L 

(r 111 (2nZ///r/'/)))iiLic 8n — 9 11 



(110) 



(r)nuc 4n2 4 

The non-relativistic limit is given by 

SC~^/C-^ = -2aZm{r)nnc ■ (HI) 

Formula (111) coincides with the related expression derived in [11] for the sphere model for 
the nuclear charge distribution, while the relativistic n-independent term in formula (110) 
differs from the corresponding term that can be derived from the formulas presented in 
[11]. Since, for the sphere model, the approach developed in [11] provides a more accu- 
rate evaluation of the nuclear size correction than the perturbation theory considered here, 
formula (110) can be improved by replacing the relativistic n-independent term with the 
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corresponding term derived from [11]. As a result, we obtain [22] 



5C-^ = ^{aZfm''{r)^A 1 + {aZf 2^(3) - V(n + 1) 



(rln(2aZmr/n))nuc 8n - 9 839 
+ 4n2 +750 



(112) 



Formulas (110) and (112) differ only by the last constant term. 

The virial relations are also helpful for calculations employing finite basis set methods or 
analytical expressions for the Coulomb-Green function. In particular, they were employed 
in [30, 31] to calculate the nuclear recoil corrections by using the B-spline method for the 
Dirac equation [32]. In that paper, using the virial relations, the original formulas for the 
recoil corrections, which contain some integrands with a singular behaviour at r — > 0, were 
expressed in terms of less singular integrals. As a result, the convergence of the numerical 
procedure for small r was significantly improved. In [8], the virial relations for diagonal 
matrix elements were used to construct Rayleigh-Schrodinger expansions for eigenvalues of 
perturbed radial Dirac equations to arbitrary order. 



7 Conclusion 

In this paper wc have considered the derivation of the virial relations for the Dirac equation 
and their applications for calculations of various physical quantities. It has been demon- 
strated that the virial relations are a very effective tool for analytical and high-precision 
numerical calciilations of the hypcrfinc splitting and the bound-electron g factor in H-like 
ions. They are also useful for calculations employing finite basis set methods and analytical 
expressions for the Coulomb-Green function. 
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